Polar phonons and intrinsic dielectric response of the ferromagnetic insulating spinel 

CdCr2S4 from first principles 
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We have studied the dielectric properties of the ferromagnetic spinel CdCr2S4 from first principles. 
Zone-center phonons and Born effective charges were calculated by frozen-phonon and Berry phase 
techniques within LSDA+U. We find that all infrared- active phonons are quite stable within the 
cubic space group. The calculated static dielectric constant agrees well with previous measurements. 
These results suggest that the recently observed anomalous dielectric behavior in CdCr2S4 is not 
due to the softening of a polar mode. We suggest further experiments to clarify this point. 

PACS numbers: 77.80.Bh, 61.50.Ks, 63.20.Dj 
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Multiferroics (MFs) displaying simultaneous ferroelec- 
tric (FE) and magnetic order are receiving considerable 
attention todayi Basic questions on the nature of ferro- 
electricity required to coexist with magnetism have been 
addressed^? facilitating the search for new materialst^ 
Much effort has been directed towards finding MFs dis- 
playing the magnetoelectric effect, due to both the desire 
to understand this particular fundamental manifestation 
of spin-lattice coupling and to the potential technological 
applications of controlling the magnetization (polariza- 
tion) by an applied electric (magnetic) fieldi^ 

One avenue for finding new MEs is to revisit known 
ferromagnetic (FM) insulators to look for ferroelectric- 
ity. CdCr2S4 is a member of a family of semiconduct- 
ing FM chalcogenide ACr2X4 spinels with A = Cd, 
Hg and X = S, Se7-^ The preponderance of evidence 
shows that CdCr2S4 crystallizes in the normal cubic 
spinel Fd3m space group over a wide temperature range, 
from 4K (lowest temperature measured) to decomposi- 
tion temperature. CdCr2S4 is a simple Heisenberg fer- 
romagnet with Cr'^+ spins (5=3/2) ordering at Tc=84K 
(0cw=155K).?^ The valence band consists of a relatively 
narrow Cr t2g peak hybridized with mostly sulfur 2p 
states, while Cr Cg states make up the lowest unoccu- 
pied states, across a gap of 1.6-1.8eV, in the conduction 
band.-W'-''^' Early interest in these materials was due to 
the observed coupling of the electronic structure and the 
lattice to the magnetic subsystemiiSiiiiiii^ While the in- 
trinsic nature of many of these effects has subsequently 
been questioned, ^^^iLiSiiS spin ordering has clearly been 
shown to have a relatively strong effect on the infrared- 
active and selected Raman-active phonon modesji^ii^ 
making the chalcogenide spinels an attractive system to 
study spin-lattice effects from first principles Recently, 
CdCr2S4 was revisited to look for ferroelectricity, lead- 
ing to the suggestion that CdCr2S4 is in fact a relaxor 
FE displaying a particularly large magnetocapacitive ef- 
fect. ^-'^■^^ A broad, frequency dependent dielectric peak 
in the real part of the dielectric constant, of the type 
associated with relaxor behavior, was observed. The di- 
electric constant at 10 Hz was shown to change by 500% 
by application of a 5 T magnetic field, with the change 
decreasing rapidly with increasing excitation frequency. 



Finally, hysteresis loops indicating nonzero polarization 
(Ps) were observed below magnetic Tc (well below the 
dielectric peak). Octahedral-site Cr'^^ off-centering was 
proposed as a possible origin of this polar behavior on the 
basis of a previous analysis^^ though such off-centering is 
unusual for magnetic ionsi^ The relaxor behavior, associ- 
ated in perovskites with polar nanoregions frustrated by 
chemical disorder,^'* was suggested to be due to geomet- 
rical frustration of local polar distortions which, as in the 
negative thermal expansion material ZrW208)2i should 
result in a very-low-frequency polar phonon mode. 

First-principles density functional methods have been 
highly successful in describing the structural, electronic, 
and magnetic properties of a variety of nonmagnetic and 
magnetic dielectrics, FEsj^i and multiferroicsiS2iS2422i2£ 
First-principles methods have emerged as an ideal tech- 
nique to differentiate between extrinsic and intrinsic 
properties of complex dielectrics and EEs, e.g. per- 
ovskite CaCu3Ti40i2.-^ In this Letter we investigate 
the intrinsic FE and dielectric behavior of cubic spinel 
CdCr2S4 from first principles. We calculate all zone- 
center phonons. Born effective charges and infrared- 
active (i.r.) mode oscillator strengths, allowing us to eval- 
uate the static dielectric response of a single-domain, sto- 
ichiometric, defect-free crystal and to evaluate the pos- 
sibility of a soft-mode-driven ferroelectric instability or 
very-low-frequency geometrically-frustrated polar mode. 

First-principles DFT calculations using PAW poten- 
tials were performed within LSDA and LSDA-I-U as 
implemented in \/ASPm^^ The wavef unctions were ex- 
panded in plane waves up to a kinetic energy cutoff of 
500 eV. Integrals over the Brillouin zone were approxi- 
mated by sums on a 6 x 6 x 6 T-centered fc-point mesh; 
for density of states calculations this grid was increased 
to 14 X 14 X 14. Phonon frequencies and eigendisplace- 
ments were calculated using the direct method where 
each symmetry adapted mode- - was moved by approx- 
imately 0.01 A. Born effective charge tensors were cal- 
culated by finite differences of the polarization using the 
modern theory of polarization**^ as implemented in VASP. 
All calculations were performed with the spin moments 
on the Cr-ions aligned ferromagnetically. 

Previous first-principles calculationsii have pointed 
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TABLE I: Structural, electronic, and magnetic properties of 
ferromagnetic CdCr2S3, Space Group: Fd3m. 



Experiment 



Theory 



structure 


300K 


40K 


LSDA 


U=3 eV 


a [A] 


10.24 


10.23 


10.24 10.06 


10.12 


u 


0.390 


0.390 


0.392 0.391 


0.390 


electronic 










E, (eV) 


1.57 


1.8 


0.77 0.00 


1.47 


magnetic 














3 


2.84 2.80 


2.88 



out the inadequacy of LSDA for investigation of the elec- 
tronic structure of CdCr2S4. In particular, LSDA fails to 
open a gap and a Cr peak lies at the edge of the va- 
lence band. This failure, and its possible implications 
for the structural instabilities and dielectric response, 
can be addressed with LSDA-|-Ui^ While the value of U 
could be calculated from first principles with constrained 
LSDA calculations, we use the alternative phenomenolog- 
ically determined value U=3 eV, which approximately re- 
produces the dominant feature of photoemission data^ 
namely a i2g peak 1.6 eV below the valence band. 

We performed a full relaxation of the lattice constant, 
a, and anion internal parameter, u, both within LSDA 
and LSDA+U. We find that the LSDA relaxed lattice 
constant, a=10.06 A, underestimates the experimental 
value by 1.8%, slightly more than that typically (<1%) 
found for LSDA calculations. Within LSDA-I-U we find a 
value a=10.12A in much better agreement (~1%) with 
the experimental value. This is consistent with other 
recent studies for LaMnOj"^ and YMnOst^S and is worth 
noting, as the LSDA+U method has only recently begun 
to be applied to structural optimization. 

The electronic structure of the fully relaxed system 
within LSDA is similar to that of previous LMTO cal- 
culations.^^ The total density-of-states (DOS) and site 
projected DOS for Cr t2g states is shown in Fig. ^ where 
a very small non-zero DOS, primarily t2g states of Cr, 
appears at the fermi level (note that this result is sensi- 
tive to volume where at the experimental lattice constant 
of a=10.24A we find a gap even with LSDA of Eg=0.77 
eV) and the t2g orbitals of the valence band occupy states 
from approximately 0.5 eV up to the Fermi level. Taking 
U = 3.0 eV opens up a gap, Eg=1.5 eV, where the top va- 
lence band is primarily oxygen p-states hybridized with 
Cr t2g (now pushed down in energy '^1.5eV) while the 
bottom of the conduction band consists of Cr e^ states 
with Cd s-states forming a wide conduction band. We 
stress, given the ad hoc procedure used to obtain U, 
these results should only be read as demonstrating that 
a reasonably good description of the electronic structure 
of CdCr2S4 is possible within LSDA-I-U. Further stud- 
ies should be performed for quantitative features of the 
electronic structure. In both LSDA and LSDA-fU, the 
calculated magnetic moments, /i=2.80/if, and ii—2.9^b re- 
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FIG. 1: LSDA (Top) Total DOS and (Middle) site-projected 
DOS on the Cr and (Bottom) eg orbitals. 
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FIG. 2: LSDA+U, U=3eV (Top) Total DOS and (Middle) 
site-projected DOS on the Cr t2g and (Bottom) Cg orbitals. 



spectively, are slightly smaller than the B^f, expected for 



C: 



.3+ 



at nominal valence . We summarized these results 



and compare with experiment in Table 

Now that we have determined that a reasonably good 
description of the structural and electronic properties of 
CdCr2S4 can be obtained within LSDA+U, we proceed 
to our main topic, the lattice dynamics and dielectric 
behavior. In Table UTI we show our calculated LSDA+U 
i.r. active TO and LO modes and compare with those 
extracted from refiectivity experiments above and below 
the FM transition. We find that all i.r. modes are sta- 
ble and in fact are quite hard, with the lowest frequency 
at 104 cm~^. Agreement of our calculated phonon fre- 
quencies with experiment is within 2%, except for mode 
T03 which still only differs by w 3.5%. To further test 
the quality of our calculations we computed the Ra- 
man active phonons. Table IIIII where again agreement 
with experiment is within a percent. This is even bet- 
ter than typical for first-principles phonon calculations, 
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TABLE II: Infrared-active phonons frequencies, lu (cm^^), 
effective plasma frequencies, fJp (cm~^), and dielectric con- 
stants of ferromagnetic CdCr2S3, Space Group: Fd3m. 



Exp. Ref. 38 Exp. Ref.J^ Theory 

300K 79K 300K 15K U=3 eV LSDA 



TO 


LO 




LO 


LO 


LO 




LO 


1 






94 


102 


104 


74 


104 


2 






239 


249 


249 


99 


247 


3 


322 


327 400 


321 


324 


339 


400 


322 


4 


377 


379 211 


378 


379 


385 


202 


366 


LO 
















1 






98 


105 


107 




106 


2 






241 


251 


251 




248 


3 


347 


350 


349 


352 


362 




349 


4 


390 


393 


392 


395 


398 




380 


^phonon 


1.9 


1.8 


2.8 


2.8 


2.3 




2.5 


Coo 


7.8±0.2 


8.0 


7.6 









where agreement to 1 or 2 wavenumbers should be con- 
sidered fortuitous. To check the possible effect of the 
underestimation of the volume on the phonon frequen- 
cies, we repeated the calculations at the experimental 
volume finding very minor changes (e.g. the lowest fre- 
quency i.r. mode is now at 102 cm~^). While non-zero 
U had a noticeable effect on the structural parameters, 
we find that the effect of non-zero U on the phonon fre- 
quencies is minimal (with an overall improved agreement 
with the experimental values) compared with the LSDA, 
also shown in Table ^ Therefore, LSDA-I-U provides a 
consistent theory to compute the lattice, electronic, and 
magnetic structures of CdCr2S4. 

The intrinsic static dielectric response is calculated as 

£0 — Eoo + / , 5— 

LO 

m ^ 

where too is the electronic dielectric constant and the 
second term on the r.h.s., tph, is the sum of contribu- 
tions from each i.r. -active phonon, called the oscillator 
strength, ^p^m and ujm are the effective plasma frequency 
and i.r. phonon frequency for mode m, respectively. A 
large dielectric response can be the result of the pres- 
ence of one or more very low frequency polar phonons 
and/or anomalously large effective plasma frequencies. 
Although can be calculated from first principles, most 
readily with density functional perturbation thcory^fl*^ 
our main interest here is the phonon contribution and 
thus we use the experimental value £00 ~ 8. We com- 
puted Vlp within LSDA-I-U from calculations of the Born 
effective charge tensors, which were found to be close 
to the nominal charges and the real-space eigendisplace- 
ments of the phonons. The results are shown in Table HII 
and compared with available experimental values, where 
again we find excellent agreement. We see that the most 
polar mode (T03), i.e. the mode with the largest 17p, 



TABLE III: Raman-active phonons frequencies, cm ^, of fer- 
romagnetic CdCr2S3, Space Group: Fd3m. 

Exp. Rcf. 18 Theory 
300K 40K LSDA U= 3 eV 



T29(l) 101 105 99 100 

Eg 256 257 257 262 

T29(2) 28 28 1 28 2 82 

T29(3) 35 1 35 3 3 5 2 352 

Aig 394 396 394 393 



is at « 300 cm^^, thus reducing its oscillator strength 
0(1). The lower frequency modes, TOl and T02, are 
weakly polar, as reflected by the very small value of Sip, 
which may be why Ref. |3^ classified these modes as non- 
polar. Also note, the relative value of fip, low for TOl 
and T02, high for T03 and T04, is due to the fact 
that the eigendisplacements of these modes are qualita- 
tively different. TOl and T02 modes consist mainly of 
Cd^"*" ions moving against both Cr'^^ and S^~ ions while 
T03 and T04 are dominated by Cr'^+ moving against 
S^~. We point out that this qualitative difference of 
the eigendisplacements is also responsible for the mode- 
dependent phonon anomalies observed at T^, which we 
will explore in a future publication.— The total contri- 
bution to the phonon part of the calculated dielectric 
constant £p/i=2.3, agrees well with that extracted from 
reflectivity measurements, eph='2,.%^ Note that £pft = l-9 
measured in Ref. 38 only contains contributions of T03 
and T04. This should be compared with our calcula- 
tions £(T03)-H£(T04)=1.7. Combining tpu with the ex- 
perimental value of £007 the computed intrinsic value of 
the static dielectric constant of CdCr2S4, £0, is found to 
be «10. This is completely consistent with the values 
obtained by reflectivity measurements and the capacitor 
measurements at 3 GHz of Ref. '2l"2?. 

Our first-principles calculations of the phonons and di- 
electric response clearly show that CdCr2S4 has neither a 
FE instability nor very-low-frequency polar modes which 
might be indicative of geometric frustration.^^ In fact, 
the intrinsic dielectric response of CdCr2S4 is found to 
be that of a normal dielectric. Further, the lack of a 
zone-center FE instability implies that no Ps is expected 
at low temperatures. These results are in direct contrast 
to previous first-principles calculations of perovskite re- 
laxors (on ordered supercells) where FE instabilities and 
anomalously large Born effective charges, yielding large 
oscillator strengths for certain modes, are found42i^ 

We considered the possibility of other structural insta- 
bilities, in particular a Cr-off-centering antiferroelcctric 
transition to the _F43m space group. In our 14-atom unit 
cell, this would be produced by the freezing in of a silent 
zone-center A211 phonon. The possibility of this phase 
transition of this type has been the subject of previous 
speculation (but never experimentally demonstrated) for 
many spinels>22ii^ Our calculation reveals that both A2U 
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modes are quite hard with frequencies greater than 300 
cm~^. In addition, the frequencies of the other silent 
zone-center phonon modes, £„, Tig, and T2u, were all 
calculated to be greater than 100 cm^^. 

Then what could be the origin of the anomalous dielec- 
tric response reported for CdCr2S4? While our results 
clearly rule out soft polar modes, one possibility might be 
the coupling of a polar mode to a zone-boundary mode, as 
we previously have suggested happens for YMnOs (note, 
unlike that of YMnOa, symmetry would require this tran- 
sition to be first order). A related mechanism might in- 
volve coupling to an incommensurate structural distor- 
tion. Although most previous experiments are inconsis- 
tent with such a picture we suggest that detailed struc- 
tural experiments should be performed. From the the- 
oretical side, first-principles calculations of the phonons 
throughout the Brillouin zone will provide further insight 
to this type of scenario. Another possibility our work sug- 
gests is that the dielectric response and polarization ob- 



served in Ref. 121123 mav not be intrinsic. While exper- 
imentally the issue of contacts have been explored, less 
has been done (to our knowledge) to characterize the de- 
fect structured- This is important given the recent claims 
that the relaxor behavior disappears in annealed single 
crystals and polycrystalline samples4i We suggest con- 
trolled defect experiments be performed to gain a more 
precise understanding of this effect. 

To summarize, the computed dielectric constant of 
10 is consistent with the measurements at frequencies 
©(GHz), while no i.r. active phonon of the cubic spinel 
structure can account for the larger responses observed 
at lower frequencies. This suggests that the origin for the 
observed relaxor behavior in CdCr2S4 is something other 
than a displacive polar soft mode. 
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